Overview
This is a tool for generating ABAQUS input data for the Concrete Damaged Plasticity (CDP) material model: compressive hardening, tension stiffening, and the compression and tension damage parameters.
Using the tool
- Compression: choose a model (Carreira & Chu, Model Code 1990, GB 50010-2010) or Test data (browse for a file, paste two columns, or load the sample compression test). A value of 0 for E, εc or ftu means “estimate it”; the value actually used is shown in blue beside each field.
- Tension: define tension stiffening as stress vs cracking strain (needs the element characteristic length lo) or stress vs cracking displacement (Linear, Bilinear or Power softening). The fracture energy Gf is entered directly or estimated from the maximum aggregate size dmax.
- Plots and tables update as you type. Each table saves as the same text file the desktop app writes, or copies as tab-separated values for Excel or the ABAQUS/CAE data table.
Compression
- Carreira & Chu (1985): σ = fcu β(ε/εc) / [β − 1 + (ε/εc)β], β = (fcu/32.4)3 + 1.55, up to ε = 0.035. fc0/fcu = 0.3.
- Model Code (1990): the CEB-FIP ascending branch up to εc,lim and the descending branch to 0.022, linear up to 0.4fcu. E = 21500((fcu + 8)/10)1/3, εc = 0.0022 and ftu from MC90 Table 2.1.2 (the ftu input is not used). Valid for fcu = 12–50 MPa.
- GB 50010-2010 (Annex C.2): σ = (1 − dc)Eε with the code's damage expressions, up to 0.022. fc0/fcu = 0.4.
- Missing E, εc and ftu are taken from EC2 Table 3.1 (except Model Code, above).
The curve is then processed for ABAQUS. The elastic limit is fc0 = (fc0/fcu) fcu, and the elastic modulus used is the secant to that point, E = fc0/εc0. Beyond it:
εin = ε − σ/E, dc = 1 − σ/fcu (after the peak; 0 before)
The hardening table stops once σ drops to 0.1fcu; the damage table starts at the peak and stops before dc reaches 0.999.
Tension
The fracture energy from dmax follows Model Code 1990: Gf = Gf0(fcu/10)0.7 with Gf0 = 0.025, 0.030, 0.058 N/mm and α = 8, 7, 5 for dmax = 8, 16, 32 mm (interpolated).
- Power (Reinhardt & Cornelissen 1984): σ/ftu = [1 + (3x)3]e−6.93x − 28xe−6.93, x = w/wc, wc = 5.14Gf/ftu. In stress–strain form εck = w/lo.
- Linear: σ = ftu(1 − w/wc), wc = 2Gf/ftu.
- Bilinear (Model Code 1990): σ = ftu(1 − 0.85w/w1) up to w1 = 2Gf/ftu − 0.15wc, then linear to zero at wc = αGf/ftu.
- GB 50010-2010: with the Chinese Code compression model, stress–strain tension stiffening uses the code's tension curve instead.
dt = 1 − σ/ftu (tables stop at σ = 0.1ftu and before dt reaches 0.95)
ABAQUS keywords
*Elastic
E, 0.2
*Concrete Damaged Plasticity
ψ, ε, fb0/fc0, K, μ (e.g. 40, 0.1, 1.16, 0.667, 0)
*Concrete Compression Hardening
<Comp. hardening: stress, inelastic strain>
*Concrete Tension Stiffening[, type=DISPLACEMENT]
<Tens. stiffening: stress, cracking strain or displacement>
*Concrete Compression Damage
<Comp. damage: dc, inelastic strain>
*Concrete Tension Damage[, type=DISPLACEMENT]
<Tens. damage: dt, cracking strain or displacement>
The plasticity parameters on the second data line are typical values for illustration (dilation angle, flow potential eccentricity, biaxial-to-uniaxial strength ratio, Kc, viscosity); this tool does not generate them.
Citation & references
Elkady, A. (2023) “ABAQUS_CDP_Generator: A tool for generating concrete damage parameters for ABAQUS”, Zenodo, Version v23.04. doi:10.5281/zenodo.7755926
- Reinhardt, H. W. and Cornelissen, H. A. W. (1984). Post-peak cyclic behaviour of concrete in uniaxial tensile and alternating tensile and compressive loading. Cement and Concrete Research, 14(2), 263–270. doi:10.1016/0008-8846(84)90113-3
- Carreira, J. D. and Chu, K.-H. (1985). Stress-strain relationship for plain concrete in compression. ACI Journal Proceedings, 82(6). doi:10.14359/10390
- CEB-FIP (1990). Model Code for Concrete Structures. Comité Euro-International du Béton and Fédération Internationale de la Précontrainte, Lausanne, 241–263.
- GB (2010). Code for Design of Concrete Structures, GB 50010-2010. National Standard of the People's Republic of China, Beijing.
- EN 1992-1-1 (2004). Eurocode 2: Design of concrete structures, Table 3.1.
Author: Dr Ahmed Elkady, University of Southampton
To report bugs or give feedback, raise an issue through the GitHub repository or e-mail a.elkady@soton.ac.uk.
Web version of ABAQUS CDP Generator v26-01. Relative to the desktop app: results update live and estimated values are shown beside the inputs instead of overwriting them; with dmax selected, Gf is always re-estimated; with test data, fcu for the Gf and EC2 estimates is the peak of the test curve.
Disclaimer
These tools are provided for research and educational purposes only, without warranty of any kind. The generated material data are estimates from the cited models and your inputs, and are only as reliable as those models within their validity ranges. Results should be checked independently by a qualified engineer before being used in design or assessment.
